Associated production of J / ^-mesons and open charm and double 
open charm production at the LHC 
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Q_i, Theoretical predictions of cross sections and properties of the J/^-meson produc- 



tion in association with an open charm hadron and formation of two open charm 
hadrons from two cc pairs in the LHC conditions are presented. Processes in both 
single and in double parton scattering mechanisms are included into consideration. 
Special attention is paid to the kinematic limits of the LHCb detector for which 
comparison with the newest experimental data is carried out. 
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b ■ I. INTRODUCTION 

In recent work of the LHCb collaboration [1] data on the double J/^-meson production 
at 7 TeV energy is presented. At first blush value of the double J /ip production cross section 
reported (5.1 ± 1.1 nb) is accordant within uncertainty limits with the predictions obtained 
in the leading order (LO) QCD calculations {2, 3|. These calculations lead to the total cross 
section value of 10 -j- 27 nb and to 3 4- 5 nb in the kinematic limits of the LHCb detector. 
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It is well known that such calculations include big uncertainties connected with the hard 
scale selection, next to the LO (NLO) contributions and allowing for the relative motion of 
c-quarks in the JM-meson. It is known that the last of this factors increases cross section 

fin 

of double quarkonia production in e + e -annihilation in several times |4H7]. 

Apart from uncertainties in the partonic cross section of double J/ip production a new 
problem arises in the LHC conditions. The hadronic cross section appears to be three orders 
of magnitude higher than the cross section of the partonic subprocess. This phenomenon 
dues to the high luminosity of low-x gluons with fraction of proton momenta of about 
10 -4 -T- 10~ 3 , which contribute most to the processes in question. Such an enhancemen; 



gave rise to discussion of double parton interactions in a single pp-collision (DPS) 18 11| 
with independent production of particles considered in each interaction. In works |12l-ll4| 
mechanism of double J/ip production in DPS approach was considered and it was shown 
that DPS can give significant contribution to the channel in question in the LHCb detector 
conditions. 

Although SPS 1 and DPS models predict somewhat different kinematic distributions for 
the J/ip pairs produced, the question if enhancement of statistics gained allows to distinguish 
this mechanisms remains open. On the other hand at least in the LO there is a qualitative 
difference between predictions for the J ~/ip + Xc production obtained in the SPS and DPS 
models. 

Moreover, additional DPS contribution should obviously express itself in other channels 
of the four c-quark domain: in the associated production of J/ip + D 2 and in the four 
.D-meson production 3 . In the beginning of 2012 first LHCb results for the channels listed 
were presented [151]. It is interesting to understand the interplay of the SPS and DPS 
mechanisms in these channels. Currently there are estimations of cross sections of the SPS 



processes contributing these final states in the LO perturbative QCD formalism [16N21]. 
In the current work we review results obtained in the LO perturbative QCD for the SPS 
contribution and estimate DPS contributions for the channels mentioned. 



1 We will address formation of considered final states in a single parton interaction as SPS. 

2 In the following we will refer to the J/ip + D production for the production of J/ip and a cc pair, from 
which at least one c-quark hadronize into an observed open charm hadron. 

3 In the following we will refer to the four D-meson production for the production of a cccc configuration, 
from which at least two c-quarks hadronize into observed open charm hadrons. 
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II. FOUR HEAVY QUARK PRODUCTION IN THE SINGLE GLUON-GLUON 

INTERACTION 

One of the first research, in which the possibility to observe four heavy quark production 
at colliders was discussed is [16] . In this work cross sections of subprocesses gg — > Q1Q1Q2Q2 
and qq — > Q1Q1Q2Q2 was estimated within LO perturbative QCD approach for the kine- 
matical conditions of the LHC and SSC. 

Slightly later the analogous processes in which quark and antiquark are bind in a doubly- 
heavy meson were investigated. Production of the S'-wave 5 c -meson in the gg — )■ B c + b + c 
and qq —> B c + b + c processes was estimated in the works |22l429|. Calculation of the P- 
wave S c -meson production cross section was done in the studies 30M32l| . These researches 



continued in the works 
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351 ] . Associated J /if) and D-meson production, as well as T and 



17H21|. 



5-meson production, was also estimated within the same technique in 

It is worth to note that doubly heavy baryon production implies production of two heavy 
quarks. Therefore, assuming that the doubly heavy baryon is created in the heavy diquark 
hadronization, one can s tudy t he doubly heavy baryon production by analogy with the 
-B c -meson production 
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36 



39]. 



Calculations show that gluon-gluon interactions provide the main contribution into 
the four heavy quark production in the LHC experiments. Quark-antiquark annihilation 
amounts to about 10%. That is why production in the gluon-gluon interactions is mainly 
discussed in this paper. 

Usually calculations are made under an assumption that initial gluons are real and their 
transverse momenta are negligible (the collinear approach). To simulate real distribution 
over the transverse momenta of initial gluons in our studies we use the Pythia 6.4 MC 



generator 



40j. In this connection it is worth noting researches [tsI . [l9 | where transverse 



momenta and virtualities of gluons are taken into account in the framework of the kr- 
factorization approach. 



III. PAIR PRODUCTION OF J/V>-MESONS IN THE LHCB DETECTOR 



Production of two charmonia in SPS can be described within perturbative QCD by the 
fourth order in as Feynman diagrams. For the J/^-mesons pair formation invariant masses 
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and quantum numbers (1 ) of two cc pairs are fixed. 

Cross section of the hard subprocess of two cc pair formation in the color-singlet (lc) 
state with m cg ~ mj/^ is proportional to 4 

where (^(O)! is the value of the cc wave function in the J/^-meson at the origin. Emergence 
of this factor dues to the approximation in which momenta of c and c quarks are parallel and 
their relative momentum is neglected in the matrix element of the subprocess in question (5- 
approximation). At large invariant masses of the J/^-meson pair cross section (TjQ) decreases 
with the rise of the full energy squared s as 

a( gg - jm*) ~ smt. (2 ) 



Numerical result of 4.1 ± 1.2 nb |2|] derived in the assumptions listed was obtained using 
the hard subprocess scale equal to the transverse mass of one of the J/^-mesons produced 
and using the CTEQ5L proton rfs fl As mentioned above this value agrees within 
uncertainty limit with the experimental value of 5.1 ± 1.1 nb measured in [1]. Variation of 
the hard scale from the one half to two transverse masses of the J/?/>-meson produced changes 
the cross section value from 5.1 nb to 3 nb. If CTEQ6LL pdfs 44j are used, cross section 
has maximum at the scale of about one transverse mass of J/tp and amounts to 3.2 nb. 
Cross sections are less at both half and double scales and are 2.8 nb and 2.6 nb respectively. 
Presence of extremum dues to the fact that with rise of the scale strong coupling constant 
decreases while gluon density grows. As well as in the manuscript [2[ we include contribution 
from the production and decay of the i[>(2S) state into the J/^-mesons yield. 

Fig. [1] shows distributions over the invariant mass of the J j ip pairs calculated within 
the assumptions mentioned for the different hard scale choices and pdf sets in comparison 
with the experimental data reported by LHCb in [l]. One can see that the shape of the 
distributions predicted is nearly the same. What concerns experimental distribution, it 
looks tilted to the bigger invariant mass values. We would like to notice that this fact 
can be accounted for by the relative c-quarks motion in the J/ip- meson. With this aim we 
calculated cross section of the process in question averaged by some "duality" region of the 



For the precise expression see [2, |4jJ, |42 1 
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m JI>pJI\In GeV 



Figure 1: Distribution over the invariant mass of the J/^-meson pairs compared with the LHCb 
measurement. Solid curves were obtained with rn^ as the hard scale, dashed — 2 • and 
dotted — 0.5 • ra^ . For every scale choice upper curve corresponds to the CTEQ5L, lower — to 
the CTEQ6LL pdf used. 



cc invariant mass: 

ff 2m D+ A d 2 a ( (cc - )f= l + (cc - )f= l) 

~ // - A f ^dm c5l dm c52 , (3) 

J J2m c dm c5l am c52 

where mo is the D-meson mass. The c-quark mass was taken equal to 

m c = 1.25 GeV. 

The A parameter can be selected in such a way that the value of the pair production cross 
section obtained coincides with the total production cross section of the J ftp J ftp, J/vptp' and 
ip'ip' final states calculated in the ^-approximation. If one takes v^S/2 for the hard scale of 
the subprocesses considered the correspondence is reached at A = 0.3 GeV: 

o d ™\gg ->■ J/^(ip')J/ip(ip'),A = 0.3 GeV) « 4.4 nb. 

At A = 0.5 GeV the cross section estimated in the duality approach is close to the value 
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Figure 2: Distribution over the invariant mass of the J/-0-meson pairs in the "duality" approach 
compared with the LHCb measurement. Solid curve was obtained with A = 0.5 GeV, dashed - 
with A = 0.3 GeV and dotted — in the ^-approximation. 

reported by the LHCb experiment: 

a dna \gg -»■ J/ip(ip')J/ip(ip'), A = 0.5 GeV) « 5.8 nb. 

Increase of A leads to the growth of the total cross section on the one hand, and improves 
agreement in the mj/^j/^ distribution on the other (Fig. [2]). 

In the LHC environment huge density of low-x gluons leads to the increase of the mul- 
tiple gluon-gluon interactions probability within one proton-proton collision. In the DPS 
approach, which implies production of particles concerned in two independent subprocesses, 
the cross section is written down as following: 



AB 
DPS 



A B 

m cr S p S cr S p g 



(4) 



2 a cS 

where the <7 e g = 14.5 mb parameter was measured in the four jets and three jets plus photon 



modes by the CDF and DO detectors 



45 



46| . The m parameter equals 1 for identical 



subprocesses and 2 for different. For the J ft/) pairs production in the LHCb conditions 
expression (j3J) leads to the following cross section value: 



a 



pp->-j/i>j/i>+x 

DPS 



4 nb. 



(5) 
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Figure 3: Distribution over the difference of the azimuthal angles of the J/^-mesons produced. 
Curve designations coincide with Fig. [TJ 

Known inclusive production cross section of the J/i/j meson in the LHCb kinematic limits, 



(Jjm = 10.5 /A), was used. In the work 14j authors note that the DPS contribution can be 
located at bigger J/ijj pair invariant masses than the SPS one. 

One of the proposed methods to distinguish the DPS signal from the SPS one is to 
study correlations between azimuthal angles of two mesons produced or between their ra- 
pidities [3, However analysis involving modelling in the Pythia generator shows that 
correlations presenting in collinear approach completely go out when including transverse 
momenta of the initial gluons into consideration (Fig. |3]). To be more precise, depending 
on the scale choice collinear or anticollinear directions of the J/ip momenta can dominate. 
Moreover, at the standard scale of one J/ip transverse mass the relative angle correlation 
is absent at all. It is the model implemented in the Pythia generator which is completely 
responsible for the distribution over the transverse momenta of the initial gluons, so model- 
independent prediction on the A0 distribution can not be derived. Investigation of the 
rapidity correlation appears to be more fruitful. In spite of narrowness of LHCb rapidity 
window (2.0 < y < 4.5) it appears to be sufficient to test QCD predictions which state 
that the difference in rapidity between the J/^-mesons produced does not exceed 2 units 
of rapidity (Fig. HJ). At the current stage DPS predicts no correlations between products of 
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Figure 4: Distribution over the difference of rapidities of the J/^-mesons produced. Curve desig- 
nations coincide with Fig. [TJ 

two partonic interactions at all. 

Apart from the correlation studies investigation of the P-wave state contributions to 
the total J/ip production can be fruitful. Indeed, section rules emerging in the CS LO 
pQCD consideration j^] imply significant limitation on the final states which can appear in 
the gluon fusion process. According to the C-parity conservation occurrence of the C-odd 
J/ipXc final state should be suppressed in the SPS. As in DPS formation of charmonia occurs 
independently, it does not have any suppression in this channel. That is why DPS should 
dominate in the J/ipXc state production (possibly followed by the xc J I ''Pi decay). It 
was pointed out in the work [13J that similar situation takes place in the J/ipT mode in 
which SPS and DPS lead to different hierarchy of the pp — > J/ip J/ip, pp — > J/if>T and 
pp — > TT cross sections. 

As it was noted, we have taken into account contribution of the pp — > J/ipip(2S) and 
pp — > ip{2S)ip{2S) processes followed by the if>(2S) — > J/ipX decays to the total J/ip pairs 
yield. ^(2S')-mesons originating from these processes can also be detected by a leptonic 
decay, just like J/ijj. It is interesting to compare ratios of J/ip J/ip and J/ipipi/lS) yields 
predicted by SPS and DPS models. 

Both in SPS and in DPS approaches ratio of different meson pair yields can be estimated 
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using values of the cc wave functions in the charmonia at the origin: 

a(pp -> J/ipJ/ip) : <r(pp ->■ J/ipip(2S)) : a(pp ->■ $(2S)i/;(2S)) « 
Vw(0) 4 : 2 ■ ^(0) 2 Vv (25) (0) 2 : </W ) 4 
« 1:1: 0.3. (6) 

where ^ J/v ,(0) = 0.21 GeV 3/2 , Vv>(2S)(°) = 0.16 GeV 3/2 . A more accurate estimate which 
allows for different meson masses leads to the relation 

a(pp -> J/ipJ/ip) : <r(pp -»■ J/ipip(2S)) : aQop -»■ tJ)(2S)iI)(2S)) « 
« 1.7:1:0.15. (7) 

Accounting contributions from the ip{2S) decays in the channels discussed one gets finally 

a(pp -> J/ipJ/tp) : cr(pp J/ipip(2S)) : <r(pp -»■ ip(2S)i/)(2S)) w 
« 2.2:1:0.13. (8) 

What concerns DPS, using inclusive J/^ and if>(2S) production cross sections equal to 



10.5 fib |48| and 1.88 fib [49] respectively, one gets 



cx(pp -> J/ipJ/ip) : J/ipip{2S)) : -»■ rp(2S)ip(2S)) 

r 2 . O ^ ~ - ~2 



2.8:1: 0.9. (9) 

It can be seen that DPS predicts slightly larger suppression of the J/iptp(2S) production 
compared to SPS. The main reason of it is that inclusive J/ip production cross section already 



includes contribution from the xc decays, which can amount up to 20 -j- 30% [50|, [5jJ. If one 
excludes expected contribution of xc decays by taking J/ip production cross section equal 
0.8 x 10.5 fib = 8.4 fib, then DPS prediction amounts to 

a(pp ->> J/ipJ/tp) : a(pp J/ipip(2S)) : a(pp -»■ ip(2S)ip(2S)) = 

= a J/i> '■ 2 ' (X J/ip Cr l(2S) '■ a i>(2S) = 

2.2:1:0.11. (10) 

Up to the uncertainties this relation coincides with the SPS prediction (JSJ). Uncertainties in 
the cross section ratios predicted by DPS can be estimated by the largest relative uncertainty 
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in the measurement of the cross sections involved. This uncertainty is maximum for the 
ip(2S) measurement and reaches 20% |49(. Unfortunately difference between relation (Q, 
which suspects feed-down from the \c production, and relations (jHJ), ( llOl) . which do not, is 
of the same order. Nonetheless it would be interesting to measure ratio of the J/ipJ/ip and 
J/i/ji/j(2S) yields experimentally. 



IV. ASSOCIATED PRODUCTION OF J/tp AND D MESON IN THE LHCB 

DETECTOR 



To compare predictions for the gg — > J/ipcc and gg — > cccc processes with experiment 
some model of the c-quark transition into a specific hadron should be used. The most 
common hardonization model is based on the assumption that charm hadron moves approx- 
imately in the same direction that the initial c quark does and obtains some fraction z of 
the quark momentum with the probability D c ^h{z) (so called fragmentation function). At 
the scale of about c-quark mass the mean z value is about 0.7. Two following parametriza- 
tions are used in our calculations: the standard parametrization of Pythia 6.4 and a pQCD 



motivated parametrization of BCFY 52j with the parameter values obtained in 53 1. 



It is worth mentioning here that as it was shown in [54H58J, there are models in which 
hadronization is not described by simple fragmentation. For example, it is reasonable to 
suppose that c-quark can pull a light quark from the sea without loosing any momentum. 
In this case it can be formally assumed that D c ^h(z) = S(z). Moreover it can be supposed 
that in some cases the final hadron momentum is even larger (by a quantity of about -^Pc) 
than the initial c-quark momentum. 

All mentioned possibilities have been considered in the present estimations of the cross 
section values. Nevertheless it should be stressed that these estimations are too rough to 
give preference to some particular hadronization model. 

Recently cross section value of the associated production of J/ip together with a D-meson 
has been measured by the LHCb collaboration for the following kinematical region jl5 |: 

• J/ip meson is produced in the rapidity region 2.0 < yj/^ < 4.0; 

• one charmed hadron is produced in rapidity region 2.0 < yo < 4.0 and has transverse 
momenta 3 GeV < < 12 GeV. 
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Calculations within the LO of pQCD lead to the cross section value 20 -=-60 nb depending 



13, 



18 



20, |21| (the scale value 



on the scale choice and the c-quark hadronization model 

was varied from y/§/4 to "*/§). Nevertheless, as it was shown is paper [27|, interaction of 
the sea c-quark from one proton with gluon from the other can essentially contribute to 
the J/^-meson and c-quark associated production, i.e. the subprocess eg — > J/ipc 5 should 
also be taken into account, as well as the main subprocess gg — > J/ipcc. It is natural 
for such an approach, that problems connected with double counting and non-zero c-quark 
mass essentially impede an accurate estimation of the calculation uncertainties. It can be 
assumed that this method is already valid at the transverse momenta of the charmed hadron 
Pt > 3 GeV ~ 2m c and that interference contributions are small. Also one can try to avoid 
double counting by subtracting the part due to the direct gluon splitting from the total 
c-quark structure function: 



fc(x,Q 2 ) = fc(x,Q 2 )- f ldZ \f X \ 2 , d X \ 2 , 2m c< Z X ) 



4:71 L Z 



- + 1-- + 



z 2 Q 2 



f g (z,Q 2 ), (ii) 



where splitting function is taken from 59]. 

The cross section value of the subprocesses eg — > J/ipc was found to be about 10 -=-40 nb. 
Therefore the contribution of such corrections to the J/ip + D associated production is of 
the same order as the contribution of the main subprocess gg — > J/ipcc. 

Thus the calculations within pQCD lead to the cross section value of about 30 -j- 100 nb 
for the J/tp + D production in the LHCb fiducial region. It should be noticed that in contrast 
to the charmonia pairs production in the associated charm production there are no C-parity 
selection rules. So one should expect not only feed-down from the J/ip + ip(2S) production 
but also from the J/ip + xc one. This contributions can increase observed J ftp + D cross 
section by up to 50%. 

In the Fig. [5] cross section distribution over the J/ip meson transverse momentum in the 
gg — > J/ip + cc subprocess is shown in comparison with the LHCb experimental data. The 
d\na/dpT distributions are plotted, i.e. spectra are normalized to unity. Both J/ijj and 
associated charmed hadron produced in the events plotted are limited to the LHCb fiducial 
region. It can be seen that at least in the high region the predicted slope is in a good 
agreement with the experimental data. It should be noticed that p'^ distribution in the 



5 Form now on summation with charge conjugate mode is implied. 
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Figure 5: Distribution over the transverse momentum of the J/^-meson in the J/ip + D production 
compared with the LHCb measurement (points for J /if) produced together with D° or D + -meson 
are shown). Solid curves were obtained with the hard scale value of 1 • m^!^ , dashed — 2 • ra^ 
and dotted — 0.5 • rrij , . Dot-dashed curve corresponds to the collinear gluon approximation. 

inclusive J / if) production measured by LHCb exhibits significantly more rapid decrease with 
the pjl^ growth. 

Cross section distribution over the D-meson transverse momentum for the same gg —> 
J/ip + cc subprocess is given in the Fig. [6] and demonstrates good agreement with the LHCb 
measurement. As in the previous figure, both spectra are normalized to unity. In contrast 
to the J/ij) signal, both predicted and measured spectra are similar to those in the inclusive 
D-meson production at LHCb |48| . 

As in double J /if? production, essential angle and rapidity correlations in the gg — > J/ipcc 
process are predicted by pQCD. Within collinear approach J/if) and D mesons move in 
the opposite directions in most cases. However no concrete prediction can be made when 
taking into account transverse gluon motion in the framework of the Pythia generator as 
the distribution is highly sensitive on the scale selection (see Fig. [7j) . 

What concerns distribution over the rapidity difference between J /if) and D-meson pro- 
duced, from Fig. [8] one can see that contrary to the two J ftp- meson production, LHCb 
rapidity window appears to be too narrow to observe rapidity correlations predicted in the 
gg — > J/ijjcc subprocess. We omit discussion of correlations between D and D mesons in the 
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Figure 6: Distribution over the transverse momentum of the D-meson in the J/ip + D produc- 
tion compared with the LHCb measurement (points for D° and D + -mesons are shown). Curve 
designations coincide with Fig. 

J/i/j + D associated production as LHCb analysis focuses on events in which one co-produced 
D meson is observed. 

The cross section value obtained by LHCb collaboration, 



a cxp (pp -> J/ip + D°(D + , D+, A+) + X) « 300 nb, 



(12) 



is several times larger than the SPS prediction of 30 -r- 100 nb. 

Let us now address to the simultaneous production of J/ip and open charm in two gluon- 
gluon interactions. Within the DPS approach cross section of the associated J/ip and D 
meson production can be expressed as follows: 



dps 

a J/tpD 



Ocff. 



(13) 



where uj/^p and <jd are cross sections of the inclusive J/ip and D-meson production in 
the LHCb acceptance correspondingly. Recalculated for the fiducial region discussed (2 < 
y < 4, > 3 GeV), these values are 9 fib and 380 fib respectively |48|, |60]. As always 
summation with the charge conjugate state is assumed. Unpublished cross section of the 
Aq inclusive production is not included in consideration. Thus the associated J/ip and D 
meson production cross section for the LHCb kinematical region within the DPS model can 



14 



1.2 - 




Figure 7: Distribution over the difference of the azimuthal angles of the J/ip and D-meson in the 
J/tp + D production. Curve designations coincide with Fig. [5j 



be estimated as 



DPS 

<J 3li>D 



240 nb. 



(14) 



As earlier, numerical value of a e g. = 14.5 mb [45j, |46( was used. One can see that DPS 
prediction is several times larger than the SPS one and within uncertainty limits agrees 
with the experimental value (TT2]) . 



V. FOUR D-MESON PRODUCTION IN THE LHCB DETECTOR 



In the same LHCb studies [15|] production of four c-quarks is investigated. Events in 
which two open charm hadrons both containing c-quark (or both containing c-quark) are 
produced in the fiducial region 2.0 < y < 4.0, 3 GeV < pt < 12 GeV were selected. 

The calculation within LO of QCD in SPS approach gives for this kinematical region 
cross section value of 

v pQCB (gg -> cccc) ~ 50 4- 500 nb 

depending on the scale selection and the c-quark hadronization model used. 

There is an indication that interaction with sea c-quarks contribute essentially into this 
process, as well as into the associated production of J / ip and c. According to our preliminary 
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Figure 8: Distribution over the difference of rapidities of the J/ip and D-meson in the J ftp + D 
production. Solid curve corresponds to the LHCb kinematic limits imposed, dashed — to the 
absence of kinematic limits. 

estimation, cross section value for the process eg — > ccc (plus charge conjugate) is about 



depending on the scale selection and the c-quark hadronization model used. The c-quark 
structure function has been taken in the form (iTTj) . 

The interactions between two sea c-quarks can also be considered. Our estimations show 
that this process can give a contribution comparable to the two processes mentioned above: 



Thus one can conclude that predictions obtained in the LO pQCD within SPS approach 
underestimate the experimental value of about 3/ib [15]. Also it is worth mentioning that 
the experimental spectra shapes also can not be exactly reproduced. 

Nevertheless some futures of the experimental spectra can be understood from such cal- 
culations using different kinematical cuts. For example the local minimum near 6 GeV in the 
experimental cross section distribution over the invariant mass m cc of two charmed particles 
is probably connected with the cut on the minimum transverse momenta at the LHCb data 



a 



pQCD (cg -> ccc) ~ 200 4- 500 nb 



a 



pQCD (cc -4 cc) ~ 40 -T- 200 nb. 
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analysis [151 ]: 



m 



loc.min 9 min 



"cc ~ ■ 

Also the rapid decrease of the cross section at m cc > 20 GeV can be explained by cut on 
the maximum transverse momenta: 



Let us now turn to the DPS contribution to the different .D-meson pairs production. 
Expression (jlj) has to be modified as experimentally observed quantities are inclusive pro- 
duction cross sections of particular types of D-mesons summed together with anti-mesons 
of the same type. These cross sections can be written down as follows: 

a? d - = a lP r + a 2 (2 P r - (p?*) 2 ), (15) 

where o\ and a 2 are cross sections of one and two cc pair production in a single proton- 
proton collision respectively and pf c is probability that c or c quark transits into detected 
hadron of type i. 

In the following we will be interested in events in which both c and c quarks form two 
.D-mesons of particular type in the detector acceptance, or it is done by pairs of identical 
quarks — cc or cc. In the first case cross section of the i type meson pair production can be 
written down as 

off* = + - (Pm 5 ) 2 + (Pf" - Ptt?/% (16) 

and in the second — as 



<r e = ^mrr + mT)(pf c - ptf) + (pf c - pit) m- (17) 



Here p\f stands for the probability for c and c quarks from one pair to transit into mesons 
of type i and j observed in the detector and probabilities for quarks from the different pairs 
are assumed independent. 

For the different i and j types of mesons analogous quantities are written down as 

<f ' = °^Tf + ^ptf - teS 5 ) 2 + iptfptf + (p?~ c - ptT - ptf) + 

+ 2pg B (pf 8 - ptf - ptf) + {pf~ c - p^ - ptfwf" - Ptf - pff)), (is) 
<* me = a 2 (0.5(ptf) 2 + 2ptfpff + 2ptf(pf 5 - ptf - ptf) + 

+ 2pt%P? 5 - Ptf ~ Ptf) + (pf 8 - Ptf ~ P-f)(Pt £ ~ Ptf ~ Ptf))- ( 19 ) 
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Mode 


diff. u 
on. ' 1 


diff- i 
&exp. : ft 


o-^ me , /ib 




D°D° 


7.2 ± 1.1 


6.2 ±0.6 


0.53 ±0.2 


0.69 ±0.07 


D°D+ 


6.0 ±0.9 


4.0 ±0.4 


0.4 ±0.1 


0.52 ±0.08 


D°Df 


2.3 ±0.4 


1.7 ±0.2 


0.16 ±0.05 


0.27 ±0.05 


D+D+ 


1.2 ±0.2 


0.78 ±0.11 


0.087 ± 0.029 


0.08 ± 0.02 


D+Df 


0.97 ±0.15 


0.55 ±0.08 


0.066 ± 0.022 


0.07 ±0.02 


DfDf 


0.19 ±0.03 




0.013 ±0.005 





Table I: Cross sections of different D-meson pairs production compared with the LHCb results. 



To solve the equations adduced we will use known inclusive production cross sections of 
particular D-meson types 48] . As LHCb collaboration presents these cross sections in bins of 
rapidity and transverse momenta, they can be recalculated into the fiducial region discussed 
(2 < y < 4, 3GeV < p^? < 12GeV). We will also assume that the total cc production cross 
section in the 7 GeV proton-proton collisions is known. It was obtained using the Pythia 
generator calibrated by known inclusive open charm production cross sections in the LHCb 



acceptance and is equal to 6.1 ± 0.9mb (48]. According to the DPS approach, cross section 



of two cc pairs production in a single proton-proton scattering is given by expression 



(T 2 



2a, 



1.3 ±0.4 mb. 



eff. 



(20) 



However until erf 1 / or aff 716 cross sections are measured there is no sufficient information 
to derive the p^ c probabilities. So we will assume that rather rigid kinematic cuts imposed 
result in the smallness of probability to observe both particles produced from a cc pair in 
the detector. Then neglecting double counting one can write down 



P 



cAc 



cVc\2 



Ptf « 2K> ) 



cVccVc 



(21) 



Under the assumptions listed equations f[T6"j) — (FUJI) can be solved. Obtained cross sec- 
tions of pair production of D°, D + and mesons are given in Table [J together with the 
values measured by the LHCb. We would like to stress one more time here that summation 
with the charge conjugate states is everywhere assumed. Generally speaking, good agree- 
ment between the DPS predictions and the experimental results is observed. Nonetheless, it 
is mentioned in Jl5| that pr-spectra of D-mesons in pair production significantly differ from 
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those in inclusive open charm production, while similar p^-behaviour could be expected in 
the DPS model. 



VI. CONCLUSION 

It is well known that the particle production multiplicity increases with the energy of 
hadronic interactions. Therefore phenomenon of multiple production should be observed for 
charmed and beauty particles as well, but at the higher energies due to the larger masses. 



At the LHC energy yield of charm particles (6.1 mb [481 ] ) is comparable to the common 
light particle yields, so production of two, three and so on pairs should be expected as well 
as single cc pair production. Recently the first data on the four c-quark production in the 
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proton-proton interactions have been obtained by the LHCb Collaboration 

From the theoretical point of view processes in single gluon-gluon interactions (such as 
gg — > cccc) are the natural source of multiple charm production. The calculations within LO 
of pQCD in SPS approach had been done earlier for the process of J/ip pair production {2} 
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421 ]. J/ip + cc associated production [16j-|2l|, and for the four c-quarks production. 



The main conclusion to be drawn from these theoretical studies and from the recent LHCb 
results is that SPS model used together with the LO pQCD can not describe all the data 
on multiple charm production. The presented analysis shows that only data on J/ip pair 
production is in satisfactory agreement with SPS LO pQCD predictions. The predictions 
obtained for the J /ip + D associated production, as well as for the four D-meson production 
underestimate the experimental data in several times. As alternative model we consider the 
simplest model of double parton scattering (DPS). In the frame work of this approach it is 
assumed that two cc pairs are produced independently in two different partonic collisions. 
DPS predictions on the cross section values fairly agree with the experimental data. As it 



was shown in 6l|, cross section of pair charm production becomes equal to the ordinary cc 
cross section at the energy of about 20 TeV. 

It is interesting to note, that for the double J/ij) production predictions of SPS and 
DPS models are fairly close, for the J/ip + D associated production the DPS prediction 
exceeds the SPS one in several times and for the four D-meson production excess is even 
higher. At first glance it seems amazing as an attempt to explain advantage of the DPS 
model by combinatorial factor only does not lead to distinction in the channels discussed. 
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Infinitesimality connected with the a s constant is same for both SPS and DPS: in SPS the 
factor is a A s and in DPS — a| x a| = otg. From our point of view the reasonable explanation 
lies in the different phase volumes for the SPS and DPS production: in SPS final state 
contains three particles for the J/ip + cc production and four for the cccc production, so 
differential cross sections of these processes peak at the larger values at the expense of 
phase volume factors. By-turn this leads to the smaller gluon luminosity as compared to 
the 2 — > 2 processes which take place in the DPS model. 

Authors would like to thank Vanya Belyaev for fruitful discussions. The work was finan- 
cially supported by Russian Foundation for Basic Research (grant #10-02-00061a), grant of 
the president of Russian Federation (^MK-3513.2012.2) and by the non-commercial foun- 
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